OPEN 3 ACCESS Freely available online 



•0-PLOS I o-^E 



CXCR7 Controls Competition for Recruitment of p- (Wt\ 
Arrestin 2 in Cells Expressing Both CXCR4 and CXCR7 cros^Mark 

Nathaniel L. Coggins^', Danielle Trakimas^^', S. Laura Chang^, Anna Ehrlich\ Paramita Ray\ 
Kathryn E. Luker\ Jennifer J. Linderman^'^*', Gary D. Luker^'^'"*' 

1 Center for Molecular Imaging, Department of Radiology, Department of Microbiology and Immunology, University of Michigan, Ann Arbor, Michigan, United States of 
America, 2 Department of Chemical Engineering, Department of Microbiology and Immunology, University of Michigan, Ann Arbor, Michigan, United States of America, 

3 Department of Biomedical Engineering, Department of Microbiology and Immunology, University of Michigan, Ann Arbor, Michigan, United States of America, 

4 Department of Microbiology and Immunology, University of Michigan, Ann Arbor, Michigan, United States of America 



Abstract 

Chemokine CXCL12 promotes growth and metastasis of more than 20 different human cancers, as well as pathogenesis of 
other common diseases. CXCL12 binds two different receptors, CXCR4 and CXCR7, both of which recruit and signal through 
the cytosolic adapter protein p-arrestin 2. Differences in CXCL12-dependent recruitment of p-arrestin 2 in cells expressing 
one or both receptors remain poorly defined. To quantitatively investigate parameters controlling association of p-arrestin 2 
with CXCR4 or CXCR7 in cells co-expressing both receptors, we used a systems biology approach combining real-time, 
multi-spectral luciferase complementation imaging with computational modeling. Cells expressing only CXCR4 maintain 
low basal association with p-arrestin 2, and CXCL12 induces a rapid, transient increase in this interaction. In contrast, cells 
expressing only CXCR7 have higher basal association with p-arrestin 2 and exhibit more gradual, prolonged recruitment of 
p-arrestin 2 in response to CXCL12. We developed and fit a data-driven computational model for association of either 
CXCR4 or CXCR7 with p-arrestin 2 in cells expressing only one type of receptor. We then experimentally validated model 
predictions that co-expression of CXCR4 and CXCR7 on the same cell substantially decreases both the magnitude and 
duration of CXCL12-regulated recruitment of p-arrestin 2 to CXCR4. Co-expression of both receptors on the same cell only 
minimally alters recruitment of p-arrestin 2 to CXCR7. In silico experiments also identified p-arrestin 2 as a limiting factor in 
cells expressing both receptors, establishing that CXCR7 wins the "competition" with CXCR4for CXCL12 and recruitment of 
p-arrestin 2. These results reveal how competition for p-arrestin 2 controls integrated responses to CXCL12 in cells 
expressing both CXCR4 and CXCR7. These results advance understanding of normal and pathologic functions of CXCL12, 
which is critical for developing effective strategies to target these pathways therapeutically. 
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Introduction 

Chemolime CXCL12 activates multiple intracellular networks, 
including mitogen activated protein kinases (MAPK), PIS kinase- 
AKT, and JAK-Stat, to control proliferation, survival, chemotaxis, 
transcription, and other cellular responses [1-3]. The numerous 
signaling pathways regulated by this chemokine correspond with 
critical functions in development, normal physiology, and disease. 
Germline deletion of CXCL12 in mice is lethal due to abnormal 
development of cardiovascular, hematopoietic, and central ner- 
vous systems [4-6] . CXCL 1 2 controls trafficking of immune cells 
and homing and retention of hematopoietic stem cells in bone 
marrow. CXCL12-dependent pathways promote growth and 
metastasis of more than 20 different human malignancies, and 
this chemokine also affects pathogenesis of other common diseases 
such as atherosclerosis, multiple sclerosis, rheumatoid arthritis and 
diabetes [7,8]. 



CXCL 12 signals through chemokine receptors CXCR4 and 
CXCR7 (recently renamed ACKR3). In cells expressing only 
CXCR4, CXCL12 binding to CXCR4 initiates signaling path- 
ways typical of seven transmembrane receptors, including 
activation of heterotrimeric G proteins and recruitment of the 
cytosolic adapter protein P-arrestin 2. The CXCR4-P-arrestin 2 
complex internalizes to endosomes, initiating P-arrestin-dependent 
signaling and ultimately leading to receptor degradation [9]. 
Conversely, CXCR7 is an atypical chemokine receptor that does 
not activate G proteins in response to CXCL12 [10]. CXCR7 
functions in part as a chemokine decoy receptor for CXCL 12, 
removing this chemokine from extracellular space and degrading it 
[11-13]. Functions of CXCR7 are enhanced by 10-fold higher 
binding affinity for CXCL 12 relative to CXCR4 and constitutive 
internalization and recycling of CXCR7 to the cell membrane 
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[12,14]. In response to CXCL12, CXCR7 also signals through P- 
arrestin 2 dependent pathways on endosomes [3,15]. 

Cells commonly co-express CXCR4 and CXCR7 under both 

normal and pathologic conditions, and studies strongly suggest 
that cells regulate levels of these receptors to respond to the 
environment and acquire new functions. For example, estrogen 
has been reported to increase expression of CXCR4 while 
reducing amounts of CXCR7 on breast cancer cells [16]. 
Activated macrophages increase mRNA and protein for CXCR7 
while downregulating CXCR4, and platelets from patients with 
acute coronary artery' disease increase CXCR7 while maintaining 
levels of CXCR4 [17,18]. In addition, tumor-initiating cells from 
some brain cancer cell lines may preferentially express CXCR4, 
contrasting with more differentiated cancer cells with greater 
expression of CXCR7 [19]. Changes in numbers of CXCR7 
versus CXCR4 receptors on cells may alter signaling pathways 
normally activated by CXCR4 alone, but reported elfects are 
contradictory [20-22]. CXCR7 has been reported to either impair 
or enhance CXCL12-CXCR4 activation of G protein signaling. 
Co-expression of CXCR4 and CXCR7 also may increase pS- 
arrestin-mediated signaling, although dynamics and distribution of 
|3-arrestin 2 between CXCR4 and CXCR7 under basal and 
Ugand-activated states remain unknown. Discordances among 
these studies with CXCR4 and CXCR7 may be due to factors 
including relative differences in ratios of CXCR4 and CXCR7 
used by different authors. 

Prior studies by our group and others have analyzed pairwise 
interactions of P-arrestin 2 with either CXCR4 or CXCR7 under 
basal conditions and in response to ligands such as CXCL12 
[14,21,23-26]. These experiments lacked the capability to 
simultaneously quantify recruitment of fi-arrestin 2 to each 
receptor in cells co-expressing both CXCR4 and CXCR7, 
precluding direct analyses of competition for this adapter protein. 
To overcome this limitation, we utilized a recendy described dual 
color click beede luciferase complementation assay for biolumi- 
nescence imaging of two different proteins interacting with a 
shared partner [27]. By fusing CXCR4 and CXCR7 to N- 
terminal fragments of click beede red and green luciferases and fi- 
arrestin 2 to the common C-terminal fragment, we could direcdy 
measure association of P-arrestin 2 with each receptor in different 
spectral windows. Dual color luciferase complementation also has 
the advantage of quantifying protein interactions in the same 
population of intact cells over time. 

To capture complex dynamics of CXCL12-dependent recruit- 
ment of P-arrestin 2 in cells expressing CXCR4 (CXCR4"*), 
CXCR7 (CXCR7'*) or both (CXCR4''-CXCR7"^, we used a 
systems biology approach combining dual color luciferase 
complementation imaging with computational modeling of 
receptor signaling and trafficking. We focused on recruitment of 
P-arrestin 2 to CXCR4 or CXCR7 since this is a common, early 
ev(-nt in ligand-dependent activation of both receptors. Based on 
imaging data from cells expressing either CXCR4 or CXCR7, we 
developed and tuned a computational model that describes 
kinetics and magnitude of interactions with P-arrestin 2. This 
approach builds on advantages of computational models to 
advance understanding of complex kinetic events and non-linear 
pathways in signaling [28,29]. Our systems biology approach 
successfully predicted perturbations in P-arrestin 2 recruitment to 
CXCR4 and CXCR7 in cells expressing both receptors and 
identified levels of P-arrestin 2 as a key control point in CXCL12- 
CXCR4 signaling. These findings elucidate how co-expression of 
CXCR4 and CXCR7 regulates P-arrestin 2 recruitment and 
underscore the power of an integrated computational modeling 



and quantitative imaging approach to investigate integrated 
functions of CXCL12, CXCR4, and CXCR7. 

Results 

Association of CXCR4 and CXCR7 with p-arrestin 2 in 
CXCR4+ and CXCR7+ cells 

To study interactions of CXCR4 and CXCR7 with P-arrestin 
2, we used a combination of data derived from a luciferase 
complementation system based on green- and red-shifted variants 
of click Ixx-tk- luciferase (Fig. lA) and outputs from a computa- 
tional model based on ordinary differential equations (Fig. IB). In 
the luciferase complementation system, CXCR4 and CXCR7 are 
fused to the N-terminal fragment of click beetie red or green 
luciferase (CBRN or CBGN), while P-arrestin 2 is fused to the 
common C-terminal enzyme fragment (CBC). We initially 
transduced MDA-MB-231 breast cancer cells with P-arrestin 2- 
CBC, so all cells express the same levels of this fusion protein. We 
then transduced cells with either CXCR4 or CXCR7 fusions. We 
validated expression of P-arrestin 2-CBC and receptor fusions by 
Western blot and RT-PCR, respectively (Fig. 2A and Table A in 
File SI). Using spectral imaging of red bioluminescence, we 
showed significantly greater basal association of P-arrestin 2-CBC 
with CXCR7-CBRN than CXCR4-CBRN (p<0.01; Fig. 2B). 

To quantify kinetics of ligand-dependent recruitment of P- 
arrestin 2 to each receptor, we treated reporter cells with 
increasing concentrations of CXCL12 from 0-1000 ng/ml and 
imaged ('(4Ls ever)' two min for 40 min and again at 90 min. We 
normalized data to \'alues from cells treated ^vith vehicle control at 
each time point. This approach focuses on rclati\'c increases in 
CXCL12-dependent recruitment of P-arrestin 2 to CXCR4 or 
CXCR7 and accounts for depletion of luciferin substrate over 
time. 

Association of CXCR4-CBRN with P-arrestin 2-CBC increased 
rapidly in a concentration-dependent manner after adding 
CXCL12 (Fig. 3 A, Fig. A in File SI). As littie as 37 ng/mL 
CXCL12 increased association of CXCR4 and P-arrestin 2 above 
basal levels, and 1000 ng/mL produced a 4-fold increase in 
bioluminescence. Interaction of CXCR4 and P-arrestin 2 peaked 
at «20-22 minutes for cells treated witii 1000 ng/mL CXCL12, 
while plateau levels occurred slightly later for cells incubated with 
lower concentrations of CXCL12. Association of CXCR4-CBRN 
with P-arrestin 2-CBC decreased in a concentration-dependent 
manner by 90 min. 

Cells expressing CXCR7-CBRN/P-arrestin 2-CBC also 
showed concentration-dependent increases in bioluminescence 
with only 4 ng/mL required to boost signal above basal levels 
(Fig. 3B, Fig. A in File SI). The difference in lower limits of 
detection for P-arrestin 2 recruitment by CXCR7 reflects greater 
affinity of this ligand-receptor pair relative to CXCL12-CXCR4 
[14]. Ligand-dependent association of CXCR7 and P-arrestin 2 
increased progressively over the fuU time course, reaching a 
maximum of 2.4-fold induction after 40 min with 1000 ng/mL 
CXCL12. Cells expressing CXCR7-CBRN exhibited more 
prolonged association with P-arrestin 2-CBC, staying constant 
through 90 min for all concentrations. After 90 min, ligand- 
dependent recruitment of P-arrestin 2 to CXCR7 was comparable 
to CXCR4. 

Describing kinetics of p-arrestin 2 recruitment and 
receptor internalization in CXCR4^ and CXCR7"^ cells by 
computational modeling 

We developed a data-driven computational model based on 
ordinary differential equations to describe receptor internalization 
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Figure 1 . Diagrams of click beetle complementation reporters and computational model of [i-arrestin 2 recruitment to CXCR4'^ or 
CXCR7*. (A) Schematic of luciferase complementation reporters for CXCR4 or CXCR7 interaction with |3-arrestin 2. (B) IVlodel schematic of receptor 
dynamics for CXCR4* cells (left) and CXCR7* cells (right) with p-arrestin 2. Note that chematic does not distinguish between endogenous p-arrestin 2 
and p-arrestin 2-CBC. 
doi:1 0.1 371 /journal.pone.0098328.g001 



and kinetics of CXCL12-dependent recruitment of P-arrestin 2 to 
either CXCR4 or CXCR7 in cells expressing a single comple- 
mentation reporter (Fig. IB, Tables 1-2; Table B in File SI). 
Consistent with experimental data, interaction of CXCR4 with P- 
arrestin 2 increases rapidly within 12-14 min, reaching a 
maximum of ~4-fold above basal levels for 1000 ng/mL 
CXCL12 (Fig. 3C). The initial increase in recruitment to CXCR4 
is due to increasing numbers of cell-surface Kgand-bound CXCR4, 
which have a higher afiinity for P-arrestin 2 than free CXCR4 
(Fig. B in File SI). More prolonged P-arrestin 2 recruitment is due 
to an increase in the number of intracellular CXCR4 molecules 
bound to P-arrestin 2. The model also shows delayed kinetics of P- 
arrestin 2 recruitment to CXCR4 with lower concentrations of 



CXCL12. Association of CXCR4 and P-arrestin 2 decreased after 
40 min for 1000 ng/mL CXCL12, declining to -75% of peak 
value by 90 min. At lower concentrations of CXCL12, recruitment 
increases at a slower rate and then plateaus. 

In contrast, association of P-arrestin 2 with CXCR7 increases 
slowly throughout 40 min for all concentrations of CXCL12 
(Fig. 3D). Interaction with P-arrestin 2 increases through 90 min in 
cells treated with 12 ng/mL or higher CXCL12, while recruit- 
ment maintains the same level by 90 min at lower concentrations. 
Similar to CXCR4, the initial increase in P-arrestin 2 recruitment 
to CXCR7 is due to an increase in the number of ligand-bound 
CXCR7 receptors on the cell surface, whereas later kinetics of P- 
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Figure 2. Luciferase complementation system reports on 
association of CXCR4 or CXCR7 witKi p-arrestin 2. (A) Expression 
of stably transduced p-arrestin 2-CBC and endogenous )3-arrestin 1/2 in 
total lysates were detected by Western blot. Blots were stripped and re- 
probed for GAPDH as a loading control. Lane 1, CXCR4-CBRN/|3-arrestin 
2-CBC; lane 2, CXCR7-CBRN/p-arrestin 2-CBC; lane 3, CXCR4-CBRN/ 
CXCR7-CBGN/|3-arrestin 2-CBC. (B) Bioluminescence in CXCR4-CBRN/P- 
arrestin 2-CBC and CXCR7-CBRN/p-arrestin 2-CBC cells was measured 
under basal conditions and 18 minutes after adding 1000 ng/ml 
CXCL12-0(. Graph shows mean values for photon flux arbitrary units 4- 
SEM for CXCR44- or CXCR74- cells (n = 4 per condition). *, significant 
difference. 

doi:1 0.1 371 /journal.pone.0098328.g002 

arrestin 2 recruitment are governed by increasing intracellular 
pools of receptors bound to P-arrestin (Fig. B in File SI). 

The model also describes experimental data for receptor 
internalization in the absence and presence of GXCL12 (Fig. 3E 
and 3F). These results establish that our computational model 
reproduces both CXCL12-dependent P-arrestin 2 recruitment to 
CXCR4 and CXCR7 and receptor internalization observed in 
dynamic biological systems. 

Predictions of p-arrestin 2 binding affinity, available p- 
arrestin 2, and available receptors 

We can use the computational model to infer mechanisms that 
drive observed behavior. In the absence of CXCL12, we calculate 
the apparent equilibrium dissociation constant of P-arrestin 2 for 
CXCR4 (Koji^ b) as over 3 times the value of the apparent 
equilibrium dissociation constant of P-arrestin 2 for CXCR7 
(.Kd,r-,,b) (1.5x10'^ nM and 4.5x10^ nM, respectively; Table 2). 
Thus, CXCR7 has higher binding to P-arrestin 2 under basal 
conditions and must recruit more P-arrestin 2 following ligand 
addition to achieve the same fold-change value as GXCR4. 

Recruitment of P-arrestin 2 to CXCR4"'" cells in response to 
1000 ng/mL CXCL12 peaks at ~ 20-22 min and plateaus or 
decreases through 90 min at all ligand concentrations. We predict 
an excess of P-arrestin 2 throughout this time (Fig. 4A). In contrast, 
the number of CXCR4 receptors able to bind P-arrestin 2 



decreases to —50% of the initial value 40 min after adding 
CXCL12 (Fig. 4B). Internahzed, ligand-bound CXCR4 is 
degraded, decreasing numbers of cell surface receptors and 
subsequently reducing the rate of P-arrestin 2 binding. Therefore, 
P-arrestin 2 recruitment in CXCR4"'' cells is limited by the number 
of ceU-surface receptors unbound to P-arrestin 2 and not the 
amount of P-arrestin 2. 

In CXCR7^ cells, CXCL12 increases association of P-arrestin 2 
with CXCR7 throughout a 40 min experiment (Fig. 4A). Our 
model predicts the number of CXCR7 receptors unbound to P- 
arrestin 2 initially decreases after adding 1000 ng/mL CXCL12 
and then partially recovers as internalized receptors recycle to the 
cell surface (Fig. 4B). Recycled CXCR7 rebinds P-arrestin 2, 
contributing to a progressive increase in interactions over time. 
Internalized CXCR7 also remains associated with P-arrestin 2, so 
complexes of CXCR7 and P-arrestin 2 accumulate intracellularly 
and continue to produce bioluminescence. 

CXCR7 decreases recruitment of p-arrestin 2 to CXCR4 in 
cells co-expressing both receptors (CXCR4^-CXCR7^) 

We next used the computational model to predict how co- 
expression of CXCR4 and CXCR7 on the same cell affects 
recruitment of P-arrestin 2 to each receptor. Simulated experi- 
ments using receptor numbers typical of our cells show that 
maximum fold change for recruitment of P-arrestin 2 to CXCR4 
decreases at all concentrations of CXCL12 in CXCR4"^-CXCR7'^ 
cells (Fig. 5A). Compared with our model output for CXCR4"'" 
cells, peak fold-induction for P-arrestin 2 recruitment to CXCR4 
decreases by -30% in CXCR4'^-CXCR7'^ ceUs with a more 
pronounced decrease over time. Conversely, the model predicts 
only a slight reduction in maximum fold-change in P-arrestin 2 
association with CXCR7 and minimal effect on progressive 
increase in signal over time (Fig. 5B). 

Our model points to a likely explanation for decreased 
recruitment of P-arrestin 2 to CXCR4 in cells co-expressing 
CXCR7. In contrast to single-receptor simulations, cell surface 
CXCR4 remains elevated through 40 min in co-expression 
simulations and does not limit P-arrestin 2 recruitment (Fig. 4B). 
Instead, the limiting factor is the amount of P-arrestin 2 available 
for binding. Free P-arrestin 2 decreases substantially in CXCR4''"- 
CXCR7^ cells, paraUehng the decrease in free P-arrestin 2 in 
CXCR7'*" cells (Fig. 4A). The decrease in free P-arrestin 2 is due to 
two factors: 1) - 50-fold greater aflTinity of CXCL12 for CXCR7 
than for CXCR4; and 2) — 8-fold greater affmity of P-arrestin 2 
for ligand-bound CXCR7 than for ligand-bound CXCR4 (see 
Table 2). These factors limit the amount of P-arrestin 2 available 
to bind CXCR4, as CXCR7 literally "steals" P-arrestin 2 away 
from the other receptor. Internalized CXCR7 remains bound to 
P-arrestin 2, further limiting amounts of free P-arrestin 2. 
Collectively, these data demonstrate that CXCR7 limits availabil- 
ity of free P-arrestin 2, thereby diminishing CXCL12-dependent 
association of CXCR4 with this scaffolding protein. 

To test model predictions, we quantified interaction of P- 
arrestin 2-CBC with CXCR7-CBGN or CXCR4-CBRN in cells 
expressing both receptors. Maximum signal for recruitment of P- 
arrestin 2 to CXCR4 increased by only —2.7-fold above control, 
representing an —35% decrease relative to cells with only CXCR4 
(Fig. 5C). Detectable recruitment of P-arrestin 2-CBC to CXCR4- 
CBRN required 111 ng/mL CXCL12, which was substantially 
greater than the concentration of 37 ng/mL CXCL12 needed to 
increase signal above baseline in CXCR4^ cells. Additionally, both 
experimental and modeling outputs showed delayed, less sustained 
recruitment of P-arrestin 2-CBC to CXCR4-CBRN in CXCR4''- 
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Figure 3. Kinetics of p-arrestin 2 recruitment to CXCR4 or CXCR7. (A and B) IVlDA-iVlB-231 breast cancer cells expressing CXCR4-CBRN/p- 
arrestin 2-CBC (A) or CXCR7-CBRN/|3-arrestin 2-CBC (B) were treated with increasing concentrations of CXCL12-0( (ng/mL) as denoted in the legend. 
Data were collected as photon flux units. Photon flux values for each time point then were normalized to values obtained for control cells not 
incubated with CXCL12 at each time point through 40 min and at 90 min. Data are expressed as mean values ± SEIVI for fold change relative to 
control (n = 4 per point). (C and D) Experimental data were used to tune parameters for a computational model describing numbers of receptors per 
cell bound to p-arrestin 2. Model outputs for CXCR4 (C) and CXCR7 (D) were plotted as fold change relative to cells not treated with CXCL12. (E, F) 
Internalization of cell surface CXCR4 (E) or CXCR7 (F) following 40 min or 30 min, respectively, of incubation with CXCL12 was measured by flow 
cytometry. Values for 0 ng/ml CXCL12 describe internalization of CXCR4 or CXCR7 in the absence of ligand. Experimental data for CXCR7 were 
replotted based on previously published results [12]. Model fits also are shown. 
doi:1 0.1 371 /journal.pone.0098328.g003 
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Table 1. Description of species included in model and steady-state values in the absence of ligand. 



Species 


Description 


Steady-state values in the absence of ligand 
for single-expressing cells* 


Steady-State values in the absence of ligand for 
co-expressing cells*" 




riee ceM-surTace L.AL.r\4 


9.5 X 1 0*^ 


1 .3 X 1 0^ 




rree ceii-surTdce t^ALK/ 


5.0 X 1 0^ 


5.1 X 1 0^ 


L,2 (nM) 


rree extracellular t^ALLiz 


0 


0 


D /M./r-nU\ 

Bg (#/cell) 


Free endogenous p-arrestin 2 


D.U X lU -n^Be Or D.UX lU — njge " r\7Bei 


c n "} D .-.1- c r\ 1 o D 
D.UX lU ~r<4Be Of D.UX lU -rijBe ~ i^JBei 




rree p-arrestin z-t^bL 


1 .5 X [SJ 


1 .5 X [B J 


"•<se (#/<:ell) 


R-4 bound to Bp 


5.8 X 1 0^ 


5.4 X 1 0^ 




R-7 bound to Bp 


6.5 X 1 0** 


6.5 X 1 0^ 


"■iBp (#/<:ell) 


R-4 bound to Bp 


8.7 X 1 0^ 


8.0 X 1 0^ 




n-7 Douno to Dp 


9.7 X 1 0** 


9.7 X 1 0^ 




R-4 bound to Lj2 


0 


0 


C7 \#/ceM) 


R-j bound to L^2 


U 


n 
U 




R-4Be bound to Lj2 


0 


0 


f~ IM./r-nU\ 


R-7Be bound to /.j2 


0 


0 


(-U/,-nU\ 

Cjsp (#/<:ell) 


R-4Bp bound to Lj2 


0 


0 


<-7sp (#/<:ell) 


R-7Bp bound to Lj2 


0 


0 


"^Be/ (#/<:ell) 


Intracellular R-4Be 


2.0 X 1 0^ 


1 .8 X 1 0^ 


n7Bei (#/<:eM) 


Intracellular R-7Be 


1 .0 X 1 0^ 


9.9 X 1 0*^ 


"•(Bp/ l#/<:eiu 


Intracellular R-4Bp 


2.9 X 1 0^ 


2.7 X 1 0^ 


"7Bp/ l#/<:eii) 


Intracellular R-jBp 


1 .5 X 1 0^ 


1 .5 X 1 0^ 


CjBei (#/<:elU 


Intracellular C^g^ 


0 


0 


Ose, (#/cell) 


Intracellular CjBe 


0 


0 


Qsp, (#/cell) 


Intracellular Qgp 


0 


0 


Osp, (#/cell) 


Intracellular CjBp 


0 


0 


R7Be„ (#/cell) 


RjBei 3fter Bg. dissociation 


2.3x10^ 


2.3x10^ 


R7Bp,7 (#/cell) 


RjBpi after Bp dissociation 


3.5x10^ 


3.5x10^ 


Qseii (#/cell) 


Qeei after Be dissociation 


0 


0 


C4Bp,7 (#/cell) 


Qsp/ after Bp dissociation 


0 


0 


Os.,, {#/cell) 


CjBei after trafficking to late 
endosomes 


0 


0 


C7Bp,7 (#/cell) 


Oep/ after trafficking to late 
endosomes 


0 


0 


L,2i (#/cell) 


Intracellular Lu 


0 


0 



•Values correspond to steady-state conditions in single-expressing cells where the total number of cell surface and intracellular receptors is 6x10 and 1.5x10 
receptors/cell for CXCR4 and CXCR7, respectively. The total number of p-arrestin 2 molecules is 5.0 x 1 0^ and 7.5 x 1 0^ molecules/cell for endogenous p-arrestIn 2 and p- 
arrestin 2-CBC, respectively. Receptor numbers are based on reasonable agreement with the data in Table C in File SI, the assumption that a large portion of the 
receptors are intracellular in the absence of ligand, and ability to fit internalization data (Fig. 3 E,F). p-arrestIn 2 numbers are based on our data suggesting that the ratio 
of probe-labeled/endogenous p-arrestin 2 is ~ 1.5 (Fig. 2A) and literature data (12). 

**Values correspond to steady-state conditions in co-expressing cells where the total number of cell surface and intracellular receptors is 6x10^ and 1.5 x 10^ receptors/ 
cell for CXCR4 and CXCR7, respectively. The total number of p-arrestin 2 molecules is 5.0 x 1 0^ and 7.5 x 1 0^ molecules/cell for endogenous p-arrestin 2 and p-arrestin 2- 
CBC, respectively. 

doi:l 0.1 371 /journal.pone.0098328.t001 



CXCR7''' cells (compare initial slopes for reporters in Fig. 3A and 
5C). 

By comparison, CXCR4 had minimal effects on interaction of 
CXCR7 with P-arrestin 2. Relative to CXCR7'^ceUs, biolumines- 
cence from CXCR7-CBGN and P-arrestin 2-CBC decreased 
minimally by 15% in CXCR4^-CXCR7"^ ceUs treated with 
1000 ng/mL CXCL12 (Fig. 5D). Co-expression of CXCR4 
modestly increased the amount of CXCL12 needed to generate 
signal for CXCR7-CBGN and P-arrestin 2-CBC from 4 ng/mL 
to 12 ng/mL. Recruitment of P-arrestin 2 to CXCR7 in dual 
reporter cells increased over 40 min, which did not differ from 
kinetics measured in cells expressing only CXCR7. 



Levels of p-arrestin 2 control CXCL12-dependent 
association with CXCR4 on CXCR4+-CXCR7^ cells 

Model predictions and experimental data show that fold-change 
in P-arrestin 2 recruitment to CXCR4 in cells co-expressing both 
receptors decreases compared with cells expressing only CXCR4. 
The model identified the amount of P-arrestin 2 as the limiting 
factor, suggesting that increasing P-arrestin 2 in CXCR4''"- 
CXCR7^ cells should alleviate suppression of P-arrestin 2 
recruitment to CXCR4. Using our model, we predicted that 
increasing P-arrestin 2-CBC by 2-fold (designated as "2x" P- 
arrestin 2) at 111 ng/mL CXCL 1 2 would sigiiificantiy prolong P- 
arrestin 2 recruitment and elevate the fold-change value (Fig. 6A), 
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Table 2. Description and values of parameters. 



Parameter 


Description 


Value 




Literature Values 


Reference 




Forward rate constant Lu binding R4 /R4Be ^f^4Bp 


2.1 xlO"^ ^' 




2.8-6.7x10"' 


[44] 


I<f.u2.7 (nM^'s"') 


Forward rate constant of Lj2 binding /R7Be /R7BP 


1.4x10"^ 




2.8-6.7x10"' 


[44-46] 


kf,B.4 ((#/cell) -h-') 


Forward rate constant of 6^ /Bp binding R4 /C4 


8.5x10"^ ■'■ (4.3x10"^ nM" 


s-T 


10"^-10"'' 


[47]*[48]" 


kf.B.7 ((#/cell) "'s"') 


Forward rate constant of or Bp binding Rj /Cj 


1.4x10"^ {7.1 xlO"^ nM" 




10"''-10"'* 


[47]*[48]" 


Kd,r4.u2 (fiM) 


Equilibrium dissociation constant of i./2 binding R4 


40 




2-27 


[49,50] 


l<D.R7.Lu (nM) 


Equilibrium dissociation constant of i.,2 binding Rj 


0.84 




0.2-0.4 


[51] 


Kd,r4b.li2 (nM) 


Equilibrium dissociation constant of Lu from /?4ee /R4Bp 


Equation (1) in text 






[38,52] 


Kd,r7b,u2 (nM) 


Equilibrium dissociation constant of /.,2 from Rjse /R7BP 


Equation (2) in text 






[38,52] 


Kd.r4.b (#/cell) 


Equilibrium dissociation constant of /Bp from R4 


7.8x10^ ^ (1.5x10^ nMr 




lO^-lO"* 


[47]'[48]» 


Kd.r7.b (#/cell) 


Equilibrium dissociation constant of /Bp from Rj 


2.3x10^ ^"^ (4.5x10^ nM)** 




10"- 10"* 


[47]"[48]'> 


Kd,c4,b (#/cell) 


Equilibrium dissociation constant of Bg /Bp from C4 


5.1 xlO^ ^ (1.0x10^ nMr 




lO^-lO^ 


[47]'[48]» 


'<D.c7,s (#/cell) 


Equilibrium dissociation constant of /Bp from Cj 


6.5x10^ "'"^ (1.3x10^ nM)'* 




10"- 10** 


[47]"[48]» 


I , -1, 
'^e,/?4S ) 


R4Be/R4Bp internalization rate constant 


2.3x10"^ 




1-2x10"' 


[50] 


I , —1, 
i^e,R7B V> 1 


RjBe /R7BP internalization rate constant 


3.9x10"^ 




1-2x10"' 


[50] 


1 1 — 1\ 
KC4B (S ) 


C4Be /C4BP internalization rate constant 


4.7x10"^ 




3x10"' 


[53] 


( 1 -1\ 
/Ce,C7e (S ) 


Oee/Oep internalization rate constant 


2.1 x10~^ 




3x10"' 


[53] 


Kff.B.4 (S"^) 


Dissociation rate constant of B^/Bp from C4Bei /C4Bpi 


7.4 xlO"'^ ' 








Kff.BJ {S ) 


Dissociation rate constant of B^Bp from RjBei ^R7Bpi 


2.5x10"^ ' ' 








. / -1\ 
^e,C7Bi V> 1 


Rate constant of trafficking of CjBei /^7Bpi to late 
endosomes 


5.5x10 " 








l^rec,R4Bi ) 


R4Bei /R4Bpi recycNng rate constant 


6.9x10"^ ■'■ 




10""-! 0"' 


[54] 


krec.R7Bii (s ^) 


RjBeii /R7Bpii recycling rate constant 


1.1 xlO"^ ■'■^ 




10""-! 0"' 


[54] 


^rec,C7Bii (s ^) 


CjBeii /^7Bpu recycling rate constant 


2.8x10"'^"'""'" 




10""-! 0"' 


[54] 


^eg,C4Bii (s ^) 


C4Beii /C4Bpii degradation rate constant 


I.OxlO"'^ 




10"^-10"'' 


[55] 


^deg,U2i (s ^) 


Lj2i degradation rate constant 


1.0 xlO"'^ *** 




10""-! 0"' 


[55] 


n4 (#/well) 


# CXCR4+ cells per well 


4.0 xlO"^ ° 








n7 {#/well} 


# CXCRZ"" cells per well 


4.0x10"^ ° 








(#/well) 


# CXCR4+-CXCR7+ cells per well 


4.0x10'* ° 










Well volume 


7.0x10"^ 









■Fit to internalization and p-arrestin 2 binding data with CXCL12 and CXCR4 in CXCR4^ cells. 
"Fit to internalization and p-arrestin 2 binding data with CXCL12 and CXCR7 in CXCR7^ cells. 

^Reference gives maximum rate of p-arrestin 2 binding as /(/^7tota/ = 0-136 s~\ Assuming a range of 10^-5x10^ p-arrestin 2 per cell gives kf~^0~^ 
Reference gives p-arrestin 2 dissociation rate constant as ^(^—0.024 s~\ which gives a p-arrestin 2/receptor equilibrium dissociation constant Kq 
**Rate constants (k) and equilibrium dissociation constants (K) are converted from #/cell to their effective value in nM using: 
1 \ ,fce/l\ /F„„ixiV,^A\ 



lO"" (#/cell) "^s" 
"lO^-lO"^ (#/cell). 



nM * s 



and 



\ L J \cellj \V,,ii^xN^y£,J 



Cell volume (Vcen) is assumed to be 8.4x10" L based on a spherical, 20 jim diameter cell. 

***These parameters do not affect model output of fold change of p-arrestin bound (see Fig. A in File SI) but are included for completeness. 

''Experimental conditions. Cells are assumed to grow to 2-3x above confluence at the time of plating. 

^Value was converted from units reported to these units by using Avogadro's number and cell volume from original paper. 

doi:l 0.1 371 /journal.pone.0098328.t002 



Cells with 2x P-arrestin 2 maintained higher levels of association 
with CXCR4 than Ix P-arrestin 2 cells through 100-min. By this 
time, cells with Ix P-arrestin 2 returned to basal levels, whereas 
cells with 2x P-arrestin 2 maintained association with CXCR4 
above baseline. 

To validate model predictions, we sorted cells for high or low 
expression of P-arrestin 2-CBG based on fluorescence from FP650. 
Western blotting showed that cells with high FP650 fluorescence 
expressed ~2-fold more P-arrestin 2-CBC than cells with low 
FP650 (Fig. 6B). Endogenous P-arrestin 2 levels were also higher in 



cells sorted for high FP650 intensity, which may be due to 
incomplete transcription of P-arrestin 2-CBC or post-translational 
cleavage of the CBC domain. However, recruitment kinetics of 
cells with low FP650 and unsorted cells showed indistinguishable 
recruitment kinetics, so parameters were conserved (Fig. C in File 
SI). Experiments with cells expressing high or low levels of P- 
arrestin 2-CBC validated overall patterns of model predictions 
(Fig. 6C; see Fig. D in File SI for raw photon flux data). For cells 
with 2x P-arrestin 2, treatment with 1 11 ng/mL CXCL12 
produced more sustained association of CXCR4 and P-arrestin 
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CXCR4 on CXCR4+ 
CXCR7 on CXCR7+ 
CXCR4 on CXCR4+-CXCR7+ 
CXCR7 on CXCR4+-CXCR7+ 



12 16 20 24 28 32 36 40 



Time (min) 



Figure 4. Modeling free p-arrestin 2 and free (unbound to p-arrestin 2) cell surface receptors over time. (A) Model output of the % of 
Initial free p-arrestin 2 through 40 min in CXCR4+, CXCR7+, or CXCR4+-CXCR7+ cells treated with 1000 ng/mL CXCL12-0(. (B) Model output of % of 
initial free (unbound to p-arrestin 2) cell surface receptors through 40 min in cells treated with 1000 ng/mL CXCL12. Legend denotes the specific 
receptor and cell type on which the receptor is expressed. 
doi:1 0.1 371 /journal.pone.0098328.g004 



2 through 100 min with significantly greater signal at 70 and 100 
min (p<0.05). Only 2x P-arrestin 2 cells also maintained 
complementation signal significantly above baseline through 100 
min. These results verify that P-arrestin 2 levels significantly affect 
both kinetics and magnitude of recruitment to CXCR4 in cells co- 
expressing CXCR7. Receptor numbers are also predicted to affect 
the kinetics and magnitude of recruitment to CXCR4 and 
CXCR7 (Fig. E in File SI), suggesting that quantitative 
modulation of P-arrestin 2 recruitment is possible. These 
conclusions hold true in an expanded sensitivity analysis (Fig. F 
in File SI) that varies CXCR4, CXCR7, and P-arrestin 2 levels 
from more physiological values (lO'^ molecules/ceU) to the lower 
limit of the overexpression system (10 molecules/cell). 

Discussion 

Precise spatial and temporal control of CXCL12 signaling is 
essential for normal development and physiology. CXCL12 
signaling regulates chemotaxis and homing of stem cells to sites 
of injury, while perturbations of CXCL12 signaling through 
CXCR4 and/ or CXCR7 drive pathogenesis of diseases such as 
cancer. Initial studies of CXCL12 signaling focused solely on 
chemokine receptor CXCR4. However, discovery of CXCR7 as a 
second receptor for CXCL12 means that biologic effects of this 
chemokine represent integrated output(s) of both CXCR4 and 
CXCR7. Prior studies indicate that functions of the CXCL12/ 
CXCR4/ CXCR7 axis are sensitive to expression of CXCR4 and 



CXCR7 on 1) separate populations of cells in the same tissue or 
organ; and 2) the same cell type [20,30,31]. Particularly for cells 
that co-express both CXCR4 and CXCR7, only limited 
information exists about how each receptor affects activation by 
CXCL12. To understand integrated functions of CXCR4 and 
CXCR7 and control these pathways for therapy, there is an unmet 
need to establish molecular mechanisms of CXCL12-dependent 
activation of one or both receptors on the same cell. 

We developed a systems biology approach to investigate 
dynamics of P-arrestin 2 recruitment to CXCR4 and/or CXCR7. 
This approach combines real-time, multi-spectral luciferase 
complementation imaging with a data-driven computational 
model. Using cells expressing complementation reporters for P- 
arrestin 2 and either CXCR4 or CXCR7, we demonstrated that 
CXCL12 caused rapid, concentration-dependent recruitment of 
P-arrestin 2 to CXCR4 that peaked within 10—20 min and slowly 
diminished through 90 min. By comparison, ligand-dependent 
interaction of CXCR7 with P-arrestin 2 increased through 90 min 
with fold-induction over basal levels comparatively less than 
CXCR4. These data are consistent with prior studies done by our 
group and others categorizing CXCR4 and CXCR7 as class A 
and B seven transmembrane receptors based on transient and 
sustained association with P-arrestin 2, respectively 
[12,14,23,32,33]. We devised and tuned model parameters using 
data from cells expressing reporters for either CXCR4 or CXCR7 
signaling. The resultant model closely reproduced the differing 
magnitude and kinetics of CXCR4 or CXCR7 association with P- 
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Figure 5. CXCR7 limits interaction of CXCR4 and p-arrestin 2 in CXCR4+-CXCR7+ ceils. (A and B) Model outputs for CXCL12-dependent 
recruitment of p-arrestin 2 specifically to CXCR4 (A) or CXCR7 (B) in CXCR4*-CXCR7* MDA-MB-231 cells. (C and D) Experimental data for recruitment of 
P-arrestin 2-CBC to CXCR4-CBRN (C) or CXCR7-CBGN (D) in CXCR4+-CXCR7+ cells. Legend shows concentrations of CXCL12-0( used for models and 
experimental data. Data were graphed as mean values ± SEIVI for fold change in bioluminescence relative to untreated cells as in Figure 1 (n = 4 per 
experimental point). 
doi:1 0.1 371 /journal.pone.0098328.g005 



arrestin 2 in cells expressing only one receptor, establishing that 
the model captures dynamics of this early step in receptor 
activation. 

Through computational modeling and experiments, we estab- 
lished that CXCR7 wins the "competition" for CXCL12- 
dependent recruitment of |3-arrestin 2 in cells that co-express 
both CXCR4 and CXCR7. Expression of CXCR7 on the same 
cells decreases the magnitude and duration of P-arrestin 2 
recruitment to CXCR4 and elevates the concentration of 
CXCL12 required to produce a signal above basal levels. By 
comparison, co-expression of CXCR4 only minimally affected 
ligand-dependent recruitment of P-arrestin 2 to CXCR7. These 
outcomes occur because CXCR7 effectively sequesters P-arrestin 2 
from CXCR4 in cells with both receptors. As predicted by 
computational modeling, increasing P-arrestin 2 partially over- 
comes suppressive effects of CXCR7 on recruitment of P-arrestin 
2 to CXCR4. These results underscore interdependent effects of 
CXCR4 and CXCR7 on responses to CXCL12 and establish P- 
arrestin 2 as a key control point in these signaling pathways [21]. 
Receptor numbers also are predicted to affect the magnitude of 
recruitment to CXCR4 and CXCR7 (Fig. E and Fig. F in File SI), 
suggesting that quantitative modulation of both the absolute 
amount and fold-change of P-arrestin 2 recruitment is possible. 



Our results for P-arrestin 2 recruitment and downstream 
signaling provide a quantitative, molecular mechanism to explain 
prior studies showing that CXCR7 may shift signaling toward P- 
arrestin 2 in cells that also express CXCR4. DecaUlot et al 
reported that co-expression of CXCR4 and CXCR7 increased co- 
immunoprecipitation of P-arrestin 2 with CXCR7, potentiating P- 
arrestin 2-dependent signaling to MAPK pathways such as 
ERKl/2 and p38 while limiting signaling mediated by G proteins 
[21]. Sierro et al also demonstrated that co-expression of CXCR4 
and CXCR7 eliminated early activation of ERKl/2 and 
produced sustained activation of these kinases in response to 
CXCL12, a characteristic feature of signaling mediated by P- 
arrestin 2 [20]. Co-expression of CXCR7 with CXCR4 also 
augmented intracellular calcium flux in response to CXCL12. Cell 
surface CXCR4 remains elevated in cells that co-express CXCR7, 
which could potentiate CXCL12-CXCR4 signaling to G proteins. 
Further modeling and experimental data are needed to establish 
effects of CXCR7 on the magnitude and duration of CXCR4 
coupling to different downstream effectors in distinct contexts. 

While our experimental and computational models include 
multiple parameters that control CXCL12 signaling, we recognize 
there are additional levels of complexity in this signaling pathway. 
We are limited by our luciferase complementation technology to 
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Figure 6. Overall levels of p-arrestin 2 limit interaction of CXCR4 and p-arrestin 2 in CXCR4'^-CXCR7'^ cells. (A to C) Figures display that 
over-expressing p-arrestin 2 increases ligand-induced p-arrestin 2-CBC recruitment to CXCR4-CBRN in CXCR4*-CXCR7* cells at extended times. (A) 
IVlodel output is plotted as fold change in the number of receptors bound to p-arrestin 2 through 100 min after treatment with 1 1 1 ng/mL CXCL12 
normalized to untreated cells at each time-point. (B) Western blot for p-arrestin 2 in cells sorted for high and low levels of fluorescence from FP650. 
GAPDH is shown as a loading control. (C) Experimental data for p-arrestin 2-CBC recruitment to CXCR4-CBRN in CXCR4*-CXCR7* cells graphed as 
mean values ± SEM for fold change of bioluminescence relative to vehicle control at 30, 70, and 100 min after treatment with 111 ng/mL CXCL12. *, 
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quantifying two pairs of protein interactions based on green and 
red spectral variants of click beede luciferase. To integrate other 
determinants of CXCL12 signaling such as chemokine or receptor 
dimers, we currently are working to incorporate additional 
complementation systems based on Gaussia or Renilla luciferases. 
These new data then will drive incorporation of additional 
elements into the computational model as needed to accurately 
describe and predict additional components of CXCL12/ 
CXCR4/CXCR7 signahng. 

Conclusions 

We have developed an integrated experimental and computa- 
tional approach to quantify, describe, predict, and validate 
dynamics of CXCL12 signaling through CXCR4 and CXCR7 
in living cells in real time. Through this approach, we have defined 
interdependent effects of CXCR4 and CXCR7 on recruitment of 
P-arrestin 2, a key node in this signal transduction pathway. In 
cells co-expressing both receptors, CXCL12 drives recruitment of 
P-arrestin 2 to CXCR7 and limits association of this scaffolding 
protein with CXCR4. Furthermore, we predicted and verified that 
amounts of P-arrestin 2 critically determine differences in 
CXCL12-association with CXCR4 versus CXCR7. Since the 
click beede luciferase complementation reporter is compatible with 
high throughput assays, these reporter cells also could be used to 
screen libraries for molecules that target rate-limiting steps in 
CXCLl 2 signaling identified by modeling. The same reporter cells 



then can be used for imaging studies in living mice, allowing us to 
refine our computational model based on in vivo data. The 
molecular imaging and mathematical systems developed in this 
work will ultimately reveal how CXCLl 2 signaling pathways 
function in normal physiology and disease and facilitate ongoing 
efforts to control these pathways therapeutically. 

Methods 

Plasmids and lentiviruses 

We used N-terminal and C-terminal fragments of click beetle 
green and red luciferases (Promega) comprising amino acids 2-4 1 3 
and 395-542, respectively, for each spectral variant [27]. We 
designated N-terminal fragments as CBGN and CBRN for click 
beetle red and green, respectively, which confer spectral charac- 
teristics of each luciferase. The common C-terminal fragment 
(CBC) complements with either N-terminal fragment. 

To sort transduced cell populations, we modified lentiviral 
vector FUGW to replace green fluorescent protein with 
mTagBFP, nuclear-localized citrine, or FP650 [34]. We cloned 
P-arrestin 2-CBC into the vector with FP650. We inserted CBGN 
fusions for CXCR4 or CXCR7 into a vector with co-expressed 
mTagBFP, and CBRN fusions were cloned into a vector with 
nuclear citrine. PGR primers used for cloning procedures are 
shown in Supplemental Methods in File SI. Amplified products 
were confirmed by DNA sequencing. 
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Cells 

We cultured MDA-MB-231 cells (ATCC) in DMEM (Life 
Technologies) with 10% serum, 1% glutamine, and 0.1% 
penicUlin/ streptomycin. We transduced 231 cells with lentiviruses 
at low multiplicity of infection for various click beede comple- 
mentation constructs as described previously [35]. We performed 
the first round of transduction with P-arrestin 2-CBC and sorted 
cells based on co-expressed FP650. Subsequent transductions 
added CXCR4 and CXCR7 fusions with CBGN, CBRN, or both. 
For cells co-expressing both receptors, we paired CXCR4-CBRN 
with CXCR7-CBGN or the reverse spectral combination. Since 
both spectral combinations performed comparably, we show data 
only for the CXCR4-CBRN and CXCR7-CBGN pair. We sorted 
transduced cell populations for mTagBFP or nuclear citrine in 
CBGN or CBRN constructs, respectively. 

Click beetle luciferase complementation for CXCR4 or 
CXCR7 interaction with p-arrestin 2 

MDA-MB-231 human breast cancer cells stably expressing 
CXCR4-CBRN, CXCR7-CBRN, CXCR4-CBRN and CXCR7- 
CBGN, or 231 control cells were seeded at 1.5x10* cells per well 
in 96 well black-wall plates. All cell lines except for 231 control 
cells also express P-arrestin 2-CBC. Cells were grown at 37°C for 2 
days before assays. We gently aspirated medium from wells and 
replaced it with 50 nL phenol red free DMEM (Life Technologies) 
with 0.2% media grade probumin (CeUiance) 30 min before 
imaging. We added 7 |a,L of a 15 mg/ mL luciferin stock and then 
incubated cells for 5 min before adding CXCL12. Immediately 
before imaging, we added 14 |iL phenol red free DMEM 
containing 0.2% probumin and increasing concentrations of 
synthetic CXCL12-a (R&D Systems). We acquired a series of 20 
images with large binning, 2 minute exposure, and open filter on 
an rVIS 100 (Perkin Ehner) for plates containing 231-CXCR4- 
CBRN, 231-CXCR7-CBRN or 231-control cells. For cells 
expressing both green and red click beetie complementation 
reporters 231-(CXCR4-CBRN)-(CXCR7-CBGN), we obtained 
20 images with large birming and 2 minute exposure, alternating 
between 530-550 nm or 690-710 nm emission filters (FVIS 200, 
Perkin Elmer). For longer time course data points, cells were 
maintained at 37°C and 5% CO2, and imaged again at 90 min. 
To determine relative induction of bioluminescence, we normal- 
ized bioluminescence for wells treated with CXCL12 to cells 
incubated with vehicle control at each time point (n = 4 per 
condition). Data were graphed as mean values ± standard error of 
the mean (SEM). 

Flow cytometry 

We analyzed cell surface CXCR4 or CXCR7 by flow cytometry 
using monoclonal antibodies 12G5 (R&D Systems) and 1 1G8 (gift 
of ChemoCentr)'x), respectively [36]. We measured receptor 
expression by mean fluorescence intensity. We performed flow 
cytometry experiments for internalization of cell surface CXCR4 
using monoclonal antibody 12G5 as described previously for 
internalization of CXCR7 [12]. Control cells were incubated 
without CXCL12 to quantify ligand-independent receptor inter- 
nalization. 

To obtain cell populations with high and low levels of P-arrestin 
2-CBC, we sorted cells by fluorescence from co-expressed FP650. 
We collected cells with the top and bottom 10% of fluorescence 
intensities. We verified that these cell populations remained stable 
by repeating flow cytometry four days later. 



Western blotting 

We analyzed endogenous P-arrestin 1 and 2 and transduced P- 
arrestin 2 in total cell lysates by Western blotting with a rabbit 
mAb (Cell Signaling) and an anti-rabbit secondary antibody 
conjugated with horse radish peroxidase (Cell Signaling). Primary 
and secondary antibody dilutions were 1:1,000 and 1:10,000. We 
detected bound antibody complexes with an ECL Plus kit 
(Amersham). 

Model 

We developed a computational model to investigate dynamics 
of P-arrestin 2 recruitment to CXCR4 and CXCR7 (Fig. IB). 
Events and pathways included are CXCL12 binding, P-arrestin 2 
recruitment, internalization, recycling, and degradation. We 
include only P-arrestin 2 because experiments show significandy 
more association between P-arrestin 2 and CXCR4 and CXCR7 
than between P-arrestin 1 and either receptor in the presence of 
CXCL12 [23]. The model includes two pools of P-arrestin 2, 
endogenous and reporter fusion to CBC, and receptors must bind 
P-arrestin 2 to internalize. We assume all internahzed CXCL12 is 
degraded and all P-arrestin 2 is recycled following dissociation 
from internalized receptors [12]. Synthesis of receptors and P- 
arrestin 2 is assumed negligible during the timescale of the 
experiments [12]. 

CXCR4 is a type A receptor, transiently binding P-arrestin 2 
[14]. In the absence of CXCL12, CXCR4 constitutively binds P- 
arrestin 2 [23]. We assume that P-arrestin 2 dissociates from 
CXCR4 not bound to ligand during internalization and that these 
receptors recycle to the cell surface. Receptors bound to ligand 
remain associated with P-arrestin 2 during internalization; 
following internahzation, P-arrestin 2 dissociates and receptors 
are routed for degradation [12,37]. 

CXCR7 is a type B receptor, tightly binding P-arrestin 2 [14]. 
Therefore, we assume that P-arrestin 2 remains associated with 
CXCR7 during internalization. AU internalized CXCR7 is 
ultimately recycled, but receptors unbound and bound to ligand 
follow distinct routes following internalization. CXCR7 not bound 
to CXCL12 is direcdy recycled. However CXCL12-bound 
CXCR7 first trafficks to late endosomes before recycling as 
CXCL12 has been shown to slow receptor recycling [12]. We 
assume that P-arrestin 2 remains bound to CXCR7 through 
trafficking to late endosomes as recycling of CXCR7 has kinetics 
similar to dissociation of P-arrestin 2 [12]. 

The mathematical model consists of coupled nonlinear ordinary 
differential equations based on mass action kinetics (Tables 1-2, 
Table B in File SI). Equations were solved using ode45 in 
MATLAB (The MadiWorks, Natick, MA). 

IVlodel Parameter Values and Initial Conditions 

We obtained initial estimates of parameter values from 
litc-raturc and our own data (Table 2). We set identical receptor 
binding and dissociation rate constants for endogenous P-arrestin 
2 and P-arrestin 2-CBC. 

Assuming no input of energy into the system, a thermodynamic 
relationship exists among the apparent equilibrium dissociation 
constants of free receptor for ligand, free receptor for P-arrestin 2, 
hgand-bound receptor for P-arrestin 2, and P-arrestin 2-bound 
receptor for ligand [38,39]: 

^^■'^^'>M2= (1) 
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D,R-jg,Li2 



Kd,R-j,Li2 ^ Kb,C-j,B 



Kd,r-j,b 



(2) 



Model simulations begin with steady-state values of aJl species in 
the absence of ligand (Table 1). We estimated total amounts of 
cell-surface CXCR4 and CXCR7 on CXCR4+, CXCR7+, and 
CXCR4'^-CXCR7"'" cells in the absence of hgand by quantitative 
receptor binding assays (Supplementary Methods and Table C in 
File SI). The number of cell surface receptors for CXCR4 and 
CXCR7 are comparable to those reported previously for cells that 
express these receptors endogenously or through gene transfer 
[40,41]. The total amount of endogenous fi-arrestin 2 was 
determined to be on the same order of magnitude as the number 
of receptors because internalization of CXCR7 proceeds through 
P-arrestin 2 binding, and experiments have shown >50% 
internalization of CXCR7 [12]. The amount of P-arrestin 2- 
CBC relative to endogenous P-arrestin was measured by Western 
blotting (Fig. 2A). 

Fitting Model to Data 

The computational model was simultaneously fit to experimen- 
tal data on both P-arrestin 2-CBC recruitment (Fig. 3A,B) and 
receptor internalization (Fig. 3E,F) resulting from: (1) CXCL12 
binding to CXCR4"^ ceUs or (2) CXCL12 binding to CXCR7"' 
cells. To compare model output to these measurements, we 
calculated the fold change in P-arrestin 2-CBC recruitment for 
CXCR4 cells at each time-point (t) and for each ligand 
concentration (L) as: 



FoldChange„,od,4{t,L) = 



R4Bpit,L)+C4Bp(t,L) + R4Bpi(t,L) 



R4Bpit,0) + R4Bpiit,0) 



(3) 



and for CXCR7"' cells as: 



FoldChangemodj(t,L) = 

RnBp it,L) + C^Bp {t,L) + Rwpiit,L) + QBpiit,L) + C-,Bpii{t,L) (4) 
RiBp(tSS) + RiBpi{tSi) 



To compare model output to experimental data on receptor 
internalization, we first initialized the model with a total number 
of receptors and P-arrestin 2, and ran the simulations in the 
absence of ligand to steady-state. This step calculated the number 
of cell surface and internalized receptors at t = 0. All cell surface 
receptors were then mathematically differentiated from internal- 
ized receptors and internalization of these cell-surface receptors 
was tracked over time. For CXCR4, this was calculated as: 



Internalization,^od.i(.t,L) = 

IU(t,L) + C4(<,L) + j?4,B„ (t.L) + R4.B, (t.L) + C4.B, (t,L) + C^_B„ (t,L) 
R,(0,L) + C4(0,L) + /?4.j), (0,L) + (0,L) + C^ji, (0,L) + Ci_Bp (0,i) 

and for CXCR7 this was calculated as: 

Internalizalion^t^jit.L) = 
R7(t,L)+ C^(t,L)+R^Js^(l,L) + R^Ji(l,L) + C^Ji,(t,L) + C^JiJt,L) 



(5) 



R^ (0,L) + CML)+ R7.Be (0,i) + (0,£) + C^,Be (0,1.) + C,^ (0,L) 



X 100% 



(6) 



where the species in these equations represent only receptors that 
are initially on the cell-surface. 

See Table 1 for symbol descriptions. 

Goodness of fit was assessed by calculating the sum of squared 
differences between model output and experimental data for P- 
arrestin 2 recruitment over time: 



SquareError - 



y^"_ J y~^™_ J {FoldChangeexp(ti,Lj) — FoldChangemod(,ti,Lj)) 



,(7) 



and between model output and experimental data for receptor 
internalization over time: 



SquareError = 



y^"_ J y^™_ J {Internalizationss_p{ti,Lj) — Internalization^odili^Lj))' 



(8) 



where n is the number of ligand concentrations tested and m is the 
number of time-points analyzed for each ligand concentration. For 
CXCR4, the square error in P-arrestin 2 recruitment was 
compared at 6 difiFerent time-points (sufficient to reproduce the 
shape of the data) between 0 and 90 minutes and for 5 different 
concentrations of CXCL12 and the square error in internalization 
was compared at 1 time-point for 5 different concentrations of 
CXCL12. This gave a total of 35 data points for model fitting to 
CXCR4 data. For CXCR7, the square error in P-arrestin 2 
recruitment was compared at 6 different time-points (sufficient to 
reproduce the shape of the data) between 0 and 90 minutes and for 
6 different concxntrations of CXCL12 and the square error in 
internalization was compared at 1 time-point for 5 different 
concentrations of CXCL12. This gave us a total of 41 points for 
model fitting to CXCR7 data. 

To find the best fit, several rounds of Latin Hypercube 
Sampling (LHS) were used to sample the parameter space [42] 
with model simulations carried out for each of 1000 parameter sets 
for each of CXCR4 and CXCR7 in each round. The total square 
error, calculated as the sum of equations (7) and (8) above, was 
calculated for each simulation. We initially varied 1 1 parameters 
for CXCR4 and 12 parameters for CXCR7, varying parameters 
+/ — an order of magnitude from the literature estimates listed in 
Table 2. The parameter set that resulted in the smallest total 
square error each for CXCR4 and CXCR7 was chosen (Table 2) 
and is used to generate the computational model portions of 
Figures 3-6. We also determined which parameters most affected 
model output using uncertainty and sensitivit)' analysis via 
calculation of Partial Rank Correlation Coefficients (PRCC) 
[28,43] (Table D in File SI). For N= 1000 runs, parameters with 
a PRCC >0.09 or <-0.09 and a p-value <0.01 were considered 
significandy different from zero. As expected, some parameters 
(e.g. binding parameters) play a major role at early time points, 
while others (e.g. internalization and recycling parameters) are 
more significant at later time points, reinforcing the need to 
include biological processes operating over multiple time frames in 
the model. 

Statistics 

We defined statistical significance as p <0.05 based on unpaired 
t-test comparisons with Welch's correction or two-way ANOVA 
on GraphPad Prism 5 software. Presented figures are represen- 
tative of at least three independent experiments for all conditions. 
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Supporting Information 

File SI Contains the following files: Figure A: Photon flux 
values for (A) CXCR4-CBRN and (B) CXCR7-CBGN corre- 
sponding to Fig 3A, 3B. Figure B: Molectjlar Species Contrib- 
uting to Beta-arrestin 2 Recruitment. Figure C: Rerrtiitment 
kinetics of 2x P-arrestin 2, Ix P-arrestin 2 and parental CXCR4^- 
CXCR7^ complementation cell lines. Figure D: Photon flux 
values for recruitment of P-arrestin 2-CBC to CXCR4-CBRN in 
2x P-arrestin 2 and Ix P-arrestin 2 cells corresponding to Fig 6C. 
Figure E: Eflect of Changing Receptor Numbers on Beta-arrestin 
Recruitment to CXCR4 and CXCR7 in Co-expressing Cells. 
Figure F: Expanded Sensitivity Analysis for P-arrestin 2 
recruitment to CXCR4 and CXCR7 in Co-expressing Cells. 
Table A: Validation of receptor expression by rjRT-PCR. Table 
B: Model Equations. Table Bl: Equations for cellular events. 
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